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We have performed numerical simulations of SU(3) gauge theory coupled to Nf — 2 flavors of 
symmetric representation fermions. The fermions are discretized with the tadpole-improved clover 
action. Our simulations are done on lattices of length L — 6, 8, and 12. In all simulation volumes 
we observe a crossover from a strongly coupled confined phase to a weak coupling deconfined phase. 
Degeneracies in screening masses, plus the behavior of the pseudoscalar decay constant, indicate 
that the deconfined phase is also a phase in which chiral symmetry is restored. The movement of 
the confinement transition as the volume is changed is consistent with avoidance of the basin of 
attraction of an infrared fixed point of the massless theory. 

PACS numbers: 11.15.Ha, 12.60.Nz, 11.20.Rd 



I. INTRODUCTION 

Many proposed extensions of the Standard Model in- 
volve new strong-coupling mechanisms that replace the 
fundamental, weakly coupled Higgs boson [1]. Among 
these are technicolor [2, 3] and "tumbling" gauge dy- 
namics [4]. The dynamics of these pictures is based on 
weak-coupling ideas: a perturbative /3 function to evolve 
the coupling constant as the energy scale falls, and most- 
attractive-channel arguments for scale separation and se- 
lection of the condensed channel in tumbling. Recently, 
groups have begun to test these proposals with nonper- 
turbative lattice methods. Most of the studies have been 
of SU(3) gauge theories with Nf > 3 fundamental fla- 
vors [5-11], but some [12-18] have considered other gauge 
groups and other representations for the fermions, where 
a richer set of phenomena might exist [19-23]. 

Some phenomenological proposals depend on novel fea- 
tures of a given gauge theory's /3 function. One possibil- 
ity is that the f3 function has an infrared-attractive fixed 
point (IRFP) at finite gauge coupling [24, 25]. The in- 
frared limit of the massless theory is then scale-invariant 
and conformal, without confinement and without spon- 
taneous breaking of chiral symmetry [26, 27]. Alterna- 
tively, the f3 function might approach zero without actu- 
ally vanishing, so that the running coupling is nearly in- 
dependent of scale over a wide range before confinement 
finally sets in at large distances. This is the scenario 
of "walking" [28, 29]. The /3 function of the massless 
theory is thus an appealing handle for nonperturbative 
study [8, 13]. 

In earlier work, we began study of the SU(3) lattice 
theory with Nf — 2 flavors of Wilson-clover fermions in 
the sextet representation [13]. Using the background field 
method (implemented in lattice work as the Schrodinger 
functional [30]), we calculated a discrete analogue of the 



/3 function on a small lattice. Contrary to perturbative 
estimates, we discovered that the discrete /3 function van- 
ishes at a fairly weak value of the renormalized coupling. 
If confirmed on larger lattices that allow extrapolation 
to the continuum limit, this would constitute an IRFP, 
indicating a conformal IR theory. The absence of sponta- 
neous chiral symmetry breaking means that this theory 
cannot be used for technicolor [31], but perhaps it can be 
incorporated in a theory of "unparticles" [32] . 

Since an IRFP implies unbroken chiral symmetry, the 
possibility of scale separation cannot arise. On the other 
hand, if the critical coupling for chiral symmetry break- 
ing is reached when the beta function is still negative, 
the fermions will condense, acquire mass dynamically, 
and decouple from lower energy scales. Only the glu- 
ons will remain, and the theory will run towards confine- 
ment. An early quenched study in the sextet theory [20] 
found large separation between the conflnement scale and 
a much higher chiral symmetry breaking scale, consis- 
tent with the larger Casimir of the sextet representation. 
When dynamical fermions are included, lattice theories 
that are similar to ours — for example, QCD with adjoint- 
representation fermions — do in fact show scale separa- 
tion, in that the critical temperature for deconfinement 
is different from (and lower than) the critical tempera- 
ture for chiral symmetry restoration [20, 22, 23]. This 
would imply, inter alia, that there is no IRFP in these 
theories. 

Even if the massless theory possesses an IRFP, how- 
ever, scale separation might reappear when one gives the 
fermions a small mass. Avoiding the IRFP, one can look 
for a remnant of a chiral transition separated from the 
confinement transition. Between the two transitions, one 
would find an intermediate phase where the phenomenol- 
ogy of broken chiral symmetry would resemble that of 
QCD even though the quarks are not confined. More 
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generally, it is interesting to see how the putative IRFP 
influences the physics of the massive theory and whether 
the latter can give evidence for or against the existence 
of the IRFP. 

We therefore present here a complementary study, 
based on straightforward calculations to determine the 
lattice theory's phase diagram and its particle spectrum 
[14]. The salient points among our results are: (1) We 
find no scale separation, meaning that there is no inter- 
mediate phase and that deconfinement is coincident with 
chiral symmetry restoration; (2) the IRFP of the massless 
theory finds support in the behavior of the deconfinement 
transition as the size of the lattice changes. 

The outline of the paper is as follows: In Sec. II we 
give a summary of our results by presenting a phase dia- 
gram of the lattice theory. This diagram shows the crit- 
ical curve Kc(/3) where the quark mass, defined through 
the Axial Ward Identity (AWI), vanishes; it shows as 
well the crossover/transition, associated with confine- 
ment and chiral symmetry, brought about by finite di- 
mensions L of the lattice. We then proceed to a presen- 
tation of our work, beginning with the definitions of the 
theory and observables in Sec. Ill and continuing to de- 
tailed results in Sec. IV. We offer further discussion in 
Sec. V. 



II. OVERVIEW 

We performed simulations on lattices of various 
sizes and measured both global observables — plaquette, 
Polyakov loop — and spectral quantities — static potential 
and hadronic correlators. If only the time dimension of 
the lattice^ Nt is finite and equal to i, then (aL)~^ can 
be thought of as a nonzero temperature; if each of the 
spatial dimensions of the lattice Ns is equal to L, then 
(aL)^ is a finite 3-volume. 



A. The deconfinement transition 

We present in Fig. 1 the phase diagram we have deter- 
mined for the lattice theory in terms of the bare gauge 
coupling P and the hopping parameter n. The solid curve 
is hic{P)i where the quark mass vanishes (see Sec. Ill 
for its definition). Kc is, in principle, a feature of the 
theory in infinite volume, though of course we have de- 
termined it on a finite lattice. 

The other (dashed) curves in Fig. 1 refer to the loca- 
tion of finite-temperature transitions, from confinement 
at small /? to non-confinement at large (3. There is one 
curve for each Nt = 6, 8, and 12. Data taken on volumes 
iV3 X Nt = 12^ X 6, 12^ X 8, and 12^ x 12 clearly show 
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^ This is distinguished from the spatial directions by the antiperi- 
odic boundary condition on the fermion field. 



FIG. 1: Phase diagram in the (/3, k) plane. The solid curve is 
Kc(/3), where vanishes; the dashed curves are Kth(/9), the 
thermal confinement transition, for three values of Nt: short 
dashes for Nt = 6, long dashes for Nt = 8 (lower curve) and 
Nt = 12 (upper curve). The star on the Kc curve marks the 
approximate location of the IR fixed point found in Ref. [13] . 
The uncertainty in «;th is in the neighborhood of 0.001. 



crossover behavior from a strong-coupling confined phase 
to a deconfined phase, as observed through the behavior 
of the Polyakov loop. The transition curves for different 
Nt appear to approach the Kc{P) curve to meet it near 
the same point. The nearest we have come to this point is 
(/3, At) = (5.1, 0.169), which is still below k = Kc- Simula- 
tion of the theory for stronger couplings is more difficult 
and we have not ventured into the /3 < 5.1 regime. 

Degeneracies among screening masses (see below) tell 
us that the deconfined phase is also a phase in which chi- 
ral symmetry is restored. We do not find two separate 
phase transitions for confinement and for chiral symme- 
try breaking. 

The finite-temperature transition might be a true 
phase transition, rounded by the less-than-infinite spa- 
tial volume, or it might be only a crossover. In QCD 
with fundamental fermions, we expect this line to be a 
line of first order phase transitions for heavy quark mass 
(small k), terminating in a second order critical point 
and becoming only a crossover as the quark mass falls. 
This is because fermions break the center Z(3) symmetry, 
and fundamental-representation fermions favor ordering 
the Polyakov loop along the positive real axis. Our data, 
however, indicate that sextet fermions break the center 
symmetry in a way that disfavors this direction and fa- 
vors argP = ±2tt/3. An example of this behavior is 
shown in Fig. 2. This means that charge conjugation 
is spontaneously broken in the high temperature phase. 
Hence, the crossover from the confined to the deconfined 
phase will become a true phase transition in the infinite 
volume limit. 

We have marked on the Kc curve the approximate lo- 
cation of the IRFP found in the massless theory, as we 
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FIG. 2: Scatter plots of real and imaginary parts of the 
Polyakov loop from simulations at f3 — 5.5, volume 8* lattices. 
Top: K = 0.140, in the confined phase. Bottom: k = 0.150, 
in the deconfined phase. 



reported in Ref. [13]. There its location was determined 
by comparing the Schrodingcr-functional effective cou- 
phng on lattices with 4^ and 8^ sites. If it is a true IRFP 
then its location will approach a limit for sufBciently large 
values of L. It is the IRFP that has caused us to focus 
our attention mainly on the region 5.1 < /? < 6.0. 

The existence of an IRFP, not its precise location, is 
what is important. The location depends on the lattice 
action and on the particular scheme used to define the 
renormalized coupling. If the theory possesses an IRFP, 
then there is a region in the space of bare coupling con- 
stants that lies in its basin of attraction — a critical region. 
Correlation functions computed at values of bare param- 
eters in this basin will show scaling (power law) behavior 
at large distance (plus, of course, cutoff-dependent be- 



havior at short distance). 

A quark mass breaks conformal invariance and so 
the critical region can include points only on the Kc(/3) 
line. Likewise, these points cannot be in the confinement 
phase: confinement provides a dynamically generated in- 
frared scale. This means that lines of finite temperature 
transitions cannot intrude into the critical region. This 
behavior is quite different from QCD with a small num- 
ber of flavors, where the deconfinement transition meets 
the Kc line at a point that moves to ever larger /? as Nt 
increases. Thus the behavior of the deconfinement lines 
is an indirect signal for or against the existence of scal- 
ing dynamics associated with an IRFP. The trend seen 
in Fig. 1 is consistent with the presence of such a critical 
region. 



B. Chiral symmetry restoration 

We have been able to distinguish two regions in the 
K < Kc plane when Nt is finite. When the coupling is 
strong {(3 small), there is confinement and spontaneous 
breaking of chiral symmetry. At weak coupling chiral 
symmetry is restored (as much as possible for Wilson 
fermions) and there is no confinement. 

A phase with chiral symmetry breaking has a pseu- 
doscalar mass m p that extrapolates to zero at zero quark 
mass as mp ~ m^. This is not the same as a strict pro- 
portionality (amp)^ (X {amq) along a line of fixed /3, be- 
cause the lattice spacing itself will change as k changes. 
Still, at small quark mass the pion mass should become 
small compared to all other dimensionful quantities, in- 
cluding the pseudoscalar decay constant fp. All these 
other quantities will be controlled by the dynamically 
generated infrared scale and will remain nonzero in the 
chiral limit. 

In a chirally restored phase we do not expect to see 
this mass hierarchy. Instead, we expect to see parity 
doubling: the pseudoscalar and scalar mesons should be- 
come degenerate, as well as the vector and axial vector 
mesons. This effect is seen in the screening masses mea- 
sured in the high-temperature phase of ordinary QCD 
[33]. In fact, we find that all four channels — scalar, pseu- 
doscalar, vector, axial vector — are close to degenerate. 
This might be interpreted as a charmonium-like spectrum 
if the quarks are heavy. If it persists as the quarks be- 
come massless then it indicates a weak quark-antiquark 
interaction, perhaps one with no boimd states. 

A naive expectation for a screening mass is that it be- 
haves as 



m 



H 



(1) 



since it/ Nt is the lowest Matsubara frequency associated 
with antiperiodic boundary conditions in a lattice of tem- 
poral length Nf We observe all of this behavior in the 
deconfined phase, and conclude that it is also a phase of 
chiral symmetry restoration. 



4 



TABLE I: Kc and mo as determined from L = 4a simulations 
with Schrodinger functional boundary conditions [13] . Linear 
interpolation may be used safely between (3 = 5.0 and 5.5 and 
between (3 = 5.5 and 6.0. 








Wo 


5.0 


.1723 


.875 


5.5 


.1654 


.887 


6.0 


.1610 


.900 


7.0 


.1536 


.916 


8.0 


.1486 


.928 



We illustrate this behavior with results from a (12 x 
8)^ X 8 volume at = 5.5, Fig. 3. The points at 
anig > 0.8 are in the confined phase; the other points 
are deconfined. Even though {amp)"^ appears to vary 
linearly with am, down to small arriq, both fp and the 
degeneracy pattern point away from spontaneous chiral 
symmetry breaking: All states become degenerate and 
afp becomes small. 

III. LATTICE ACTION AND SIMULATION 
DETAILS 

A. Lattice action 

Our lattice theory is defined by the single-plaquette 
gauge action and a Wilson fermion action with added 
clover term [34]. We modify the clover term's coefficient 
via tadpole improvement by setting csw — ^/uq. We 
adopt the conventional self- consistent determination of 
the tadpole improvement factor, viz., 

4 = l{TrUp)^^, (2) 

where Up is the usual plaquettc variable. This is 
the same action that we used in Ref. [13]. In that 
paper we determined Kc on lattices with volume 4^, 
with Schrodinger- functional boundary conditions, and we 
fixed Uq according to the space-space plaquette averages 
on those lattices at k = Kc (see Table I). Rather than 
recalculate uq at every (/?, k) , we choose to define the the- 
ory at every /3 by fixing uq to take the value in Table I, 
irrespective of k. 

B. Data sets 

All simulations used the standard hybrid Monte Carlo 
algorithm. The trajectories in various runs were of 
lengths between 0.5 and 1.0, and the time steps ranged 
from 0.02 (at heavy quark masses) to 0.005 (for light 
masses). The data sets at each of our (/?, k) values con- 
sist of 300 to 1000 trajectories, with every fifth trajectory 
used for spectroscopy. 

Results from five different simulation volumes are re- 
ported in this work: 



• 12^ X 6 and 12^ X 8 are conventional finite- 
temperature lattices. We use them to find Kth(/3) 
curves. 

• (12 X 8^) X 8 allowed faster runs than 12'^ x 8 
and showed the same finite-temperature physics, 
though transitions are rounded by the smaller spa- 
tial volume. 

• 8^ X 12 is a "zero temperature" lattice compared 
to iVt = 8. We use it to study how the spatial size 
L = 8 intrudes on the qq potential and on meson 
masses. 

• 12^ has two roles. One is as a "zero temperature" 

lattice, as long as we stay on the strong-coupling 
side of the TVf = 12 confinement transition. The 
other role is as a finite-temperature lattice that per- 
mits us to observe directly the movement of the 
Kc(/3) curve when Nt changes from 6 to 8 to 12. 
(We are well aware that systematic thermodynam- 
ics studies require A^^ x Nt lattices with Ns ^ Nt, 
so we will be careful not to over-analyze our re- 
sults.) 

We also performed exploratory runs with 8^ volumes. 

Dimensions of size 12 are where we determine meson 
masses. If the dimension is temporal then the masses 
are conventional spectroscopic masses; if spatial then the 
masses are screening masses, affected by the Matsubara 
frequencies that create non-zero momentum transverse 
to the meson propagation. Lattices where Nt = 12 are 
also where we calculate the (spatial) qq potential from 
Wilson loops. 

A map showing the bare parameters of our runs is 
shown in Fig. 4. We concentrated on two areas of the 
(/3, k) plane: The upper band of points covers the area 
near Kc, while the lower band of points is in the vicinity 
of the deconfinement transitions. 



C. Meson propagators 

In order to measure masses when Nt = 12, we effec- 
tively double Nf with a trick that has been used by sev- 
eral groups for computing weak matrix elements [35-38] . 
This combines periodic with antiperiodic boundary con- 
ditions as follows. Take a valence Dirac operator with 
periodic temporal boundary conditions and compute its 
propagator, ^^(a;) (we assume a source at t = for sim- 
plicity). Take a second valence Dirac operator with an- 
tiperiodic temporal boundary conditions, and compute 
its propagator Sa{x). Now add the propagators to pro- 
duce 
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FIG. 3: Screening masses and fp for /3 = 5.5 on volume (12 x 8^) x 8. In (a) wc plot the squares of the quantities, while in 
(b) we plot the quantities themselves. Crosses show fp, pseudoscalars are diamonds, vectors are squares, octagons are axial 
vectors and bursts are scalars. 
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FIG. 4: Map of our runs in the (/?, k) plane. The upper set 
of points parallels the Kc line in the deconfined phase. The 
lower set of points populates the regions of Kth lines. 



and use this propagator to construct meson correlators, 
e-g-, 



(4) 



for a pseudoscalar meson. The resulting correlator will 

be a hyperbolic cosine with midpoint at i = 7V(. (We use 
the same method in measuring screening masses when 



= 12.) 



D. Meson masses, decay constant, and quark mass 

After this, our fitting methodology is standard. We 
perform spectroscopy using quark propagators computed 
in Coulomb gauge, generated from Gaussian wall sources. 
Wc take masses from correlated fits to long temporal 
ranges and use a combination of stable effective mass 
fits, high confidence levels, and stability of our masses 
under variation of fit range to choose a result and its 
uncertainty. As a general rule, when a meson is heavy, 
single-exponential fits to correlation functions work well, 
but as the meson mass falls below auriH ^ 1 we need 
two-exponential fits to give a stable lightest mass. 

We also measured the pseudoscalar decay constant and 
the AWI quark mass. The latter is defined through 

dtY,{Mx,t)Xm =2m,Y,{P{x,t)X{Q)) . (5) 

X X 

where Aq = ip"/o"/5ip and P = V'75'0- For consistency 
with the conventions of our Schrodinger functional cal- 
culation, the derivative is taken to be the naive differ- 
ence operator 9^/(a;) = [f{x + fi) — f{x — (We set 
a = 1 in this discussion.) The source ^(0) is arbitrary; 
for making the measurement as part of a spectroscopic 
study we use one made of our Coulomb gauge Gaussian 
wave functions. When it is chosen to have pseudoscalar 
quantum numbers, and when we assume that a single ex- 
ponential dominates the correlator, each side of Eq. (5) 
will be proportional to cosh(TO7r(t — 2iVt)). We then ex- 
tract the quark mass by performing a three-parameter 
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{Z, mp, ruq) correlated fit to 

X 

and 

y {Aoix,t)X{0)) = - .^f""^ fe-"^* - e-™^(2iv,-t)^ _ 
smnrnp V / 

(7) 

If we were to convert the lattice-regulated AWI quark 
mass to some continuum regularization, it would require 
multiplicative renormalization. We will neglect this over- 
all renormalization factor in our discussion. For Wil- 
son fermions, radiative effects mix the local axial current 
= i'lfilb'^ with 9^ ■^75'!/'. The mixing is small in per- 
turbation theory [39], so wc neglect it as well. 

We determine the pseudoscalar decay constant from 
the matrix element of the axial vector current, mpfp = 
(0|j4o|P). Specifically, wc perform a thrcc-paramctcr cor- 
related fit to the correlator of a Gaussian source and sink, 

{X{x,t)X{0)) = (e-™^* + e-™^(2jv,-t)j ^ 

X 

(8) 

and to the correlator of a Gaussian source and an axial 
current sink, 

J2 {Ao{x,t)X{0)) = ^ (e-™^* - e-™^(2jv.-t)j _ 

X 

A lattice determination of a continuum decay constant 
is a little involved. Because we are using a non-chiral lat- 
tice action and a matrix element that does not precisely 
realize a Ward identity, there is a lattice-to-continuum 
conversion factor Za between the lattice matrix element 
and a continuum-regulated decay constant, 

fr* = zaIp. (10) 

We are aware of two ways to compute Za- One is nonper- 
turbative, through the RI (regularization independent) 
scheme. The other is through perturbation theory. For 
our exploratory study, we believe perturbation theory is 
adequate. In the context of tadpole-improved perturba- 
tion theory [40], 

C^iR) is the appropriate Casimir (4/3 for fundamen- 
tals, 10/3 for sextets), W is & numerical factor (5.79 
for our tadpole- improved clover action), and the "tad- 
pole factor" ^1 — corrects the field renormaliza- 
tion of lattice Wilson fermions compared to the con- 
tinuum case. Tadpole-improved perturbation theory is 
designed to pick a scale for the coupling (the so- 
called q* scale in its jargon). Again for our action. 



q* = 2.41/a. Computing the coupling from the plaque- 
tte (g2(g = 3.41/a) = -3 ln(Tr J7p/3)) and running to 
the desired q*, we can estimate Za- The bottom line 
is that the non-tadpole-improvement part of Za (every- 
thing except the (1 — ^) factor) is essentially constant, 
equal to about 0.8 over our parameter set. For a quantity 
which is supposed to be unity plus perturbative correc- 
tions, the deviation is uncomfortably large. If we were 
going to make phenomcnological use of f-^ - for example, 
to predict the technirho mass from a lattice measurement 
of mp/ f.^ in a QCD-like theory — this would be cause for 
concern. Since our only interest here is in whether fp 
vanishes in the chiral limit, this is probably unimportant. 



E. The Kc curve 

The Kc curve shown in Fig. 1 (see also Table I) was 
determined in [13] by demanding rriq = on a lattice 
with 4^ sites and Schrodinger functional boundary con- 
ditions. The matrix elements in the AWI (5) were evalu- 
ated ait = 2 with the operator X{Q) defined through the 
boundary conditions themselves. Finite- volume effects 
can shift mq{K) and thus Kc{I3). We check this by calcu- 
lating mq{K) via the spectroscopic fit on (12 x 8^) x 8 and 
8^ X 12 lattices (Fig. 5). Apart from the data at strongest 
coupling {(3 = 5.1), the curves are quite smooth, and on 
extrapolating to niq = we find a shift upward of less 
than 2 x 10^^ in Kc, roughly the thickness of the solid 
curve in Fig. 1. 



F. Heavy quark potential 

In measuring the static potential between fundamen- 
tal representation sources, our approach is, again, fairly 
standard. We extract V{r) from the effective masses 
of i? X T Wilson loops after one level of HYP smear- 
ing [41, 42]. The short-distance effects of the HYP smear- 
ing are corrected using a fit to the perturbative lattice 
artifacts. With fit parameters A, B, a, and a, we write 
V{r) = Kont + aAViatt, where 

Vcont = - + B + ar (12) 
r 

and 

= X^pert(r-) - I - (13) 

A^att is the difference between the exact lowest-order 
perturbative lattice propagator and 1/r. We quote val- 
ues for the string tension a^a and the Sommer parameter 
ro/a [43] . For successful determination of the potential 
we demand consistent fit parameters from several differ- 
ent temporal sizes of the Wilson loops — fits toT = t and 
t + 1 for t = 3, 4:, and perhaps 5, for example. 
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FIG. 5: AWI quark mass from (12 x 8^) x 8 lattices. Curves 
are for (right to left) (3 = 5.1, 5.2, 5.5, 5.7, 6.0, and 7.0. 



IV. SIMULATION RESULTS 
A. Identifying the deconfinement line 

The straightforward way to find where the confine- 
ment scale crosses some distance L is to find the finite- 
temperature confinement transition/crossover curve for 
= L. This wc do via a conventional Polyakov loop 
analysis, looking for metastability, phase coexistence, or 
a rapid variation in (P). Fig. 1 shows the locations of 
the confinement transitions for the three temporal sizes. 
Nt ~ 6, 8, and 12. We have data ranging from heavy to 
light quark masses. As we approach the Kc curve or make 
/3 smaller than about 5.1, the number of conjugate gra- 
dient steps needed to compute the fermionic force grows 
steeply. Hence we are unable to say whether the transi- 
tions actually join the Kc line. 

Measurements of the static potential on lattices with 
large Nt agree well with Polyakov loop data. In the con- 
fining phase, we use a confinement scale Rc such as the 
Sommer radius ro or l/\/a- (for string tension a) to define 
a length scale. Then lattices with a spatial size L that 
is large compared to Rc exhibit confinement, and those 
with L < Rc do not. The ratio Rc/a tends to grow as 
gauge coupling P increases and as k rises. 

Examples of two fits to the potential are shown in Fig. 
6. On the left, we show results from an 8^ x 12 lattice 
at /3 = 5.3, K = 0.155. The different plotting symbols 
denote two different temporal sizes of the Wilson loops. 
This potential is clearly confining, with a large value of 
a^a and small Sommer parameter ro/a, either of which 
can be used to define Rc/a. Moving to larger /3 and/or k 
will enlarge Rc until confinement will no longer be evident 
for L = 8. 



The right panel is from a 12"* lattice at /3 = 5.7, 
K = 0.158. The potential has clearly flattened: a^a has 
become small and ro/a has grown. The different plotting 
symbols show that our results are no longer independent 
of the temporal size of the Wilson loop, and we are un- 
able to quote any number for a^cr or ro/a. Here, then, 
even L = 12 cannot contain the confining region of the 
potential. 

At very small (3 we move deep into the strong coupling 
phase. The Sommer parameter ro/a becomes small and 
the fitter can no longer determine it. The potential is 
purely linear, V{r) ^ err. This means that, for all practi- 
cal purposes, there is a rather narrow region of bare pa- 
rameters for which we can determine a^a and ro/a from 
our simulation volumes. This region coincides roughly 
with the vicinity of the Nt = 8 deconfinement transition. 



B. Hadronic spectroscopy 

As we have described, we use the behavior of fp and 
the P, V, A, and S masses to characterize the chiral 
symmetry aspects of each phase. If a decrease in the 
quark mass induces a drop in mp relative to the other 
masses while fp shows little variation, we call the phase 
chirally broken.^ If on the other hand we observe parity 
doubling — degeneracy of the P and S masses, and of the 
V and A masses — along with an fp that falls toward 
zero along with arUq, we conclude that chiral symmetry 
is restored. 

We expect the strong coupling phase to be a phase in 
which chiral symmetry is broken [44]. We find that in 
this phase amp, amy, and afp are easy to extract. As 
most of our data are taken at fairly large quark masses, 
it is hard to say more than that amy > amp. The P 
and V masses fall as k grows. The S and A signals are 
poor. Their masses are large — greater than amy — and 
the fits are unstable. 

The P mass is not affected much by the deconfinement 
transition (apart from /3 = 5.1). Along most of the tran- 
sition line, particle masses are large, mna ^ 1. 

The P and V masses continue to fall as k is taken closer 
to Kc (i.e., as amq is taken to zero). As we move into the 
deconfined phase, the signals in the S and A channels im- 
prove and their masses fall as well. All four masses — P, 
S, V, and A — become nearly degenerate. At the same 
time, afp becomes much smaller than any of the meson 
masses. As we mentioned above, the simplest interpre- 
tation of our data is that chiral symmetry is restored in 
the deconfined phase. The fourfold degeneracy may be 
reminiscent of the quark model with a weak hyperfine 
interaction, but another way to look at it is as a world 



^ For confirmation, wc could then ask if (mp/A)^ oc niq/A for 
some choice of scale-setting mass A; but to do this successfully 
amounts to finding a curve of constant physics. See below. 
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FIG. 6; Examples of potentials from our data sets. On the left, the string tension is large, ro/a is not too small, and we can 
obtain a good fit to V{r). On the right, the string tension has become small and we cannot perform a reliable fit to V{r). Left 
panel, 8^ x 12, /? = 5.3, k = 0.155. Right panel, 12^, /3 = 5.7, k = 0.158. 



where quarks form no bound states at all. This, too, is 
characteristic of a massless limit with exact, unbroken 
chiral symmetry [44]. 

We illustrate this result with several sets of spectra 
(Figs. 3 and 7-11). In all cases the right panel displays a 
massive quantity (in lattice units, orn/f) versus the AWI 
quark mass am^, while the left panel displays the squared 
quantity versus aniq. In all graphs, crosses show the de- 
cay constant (with tadpole factor), while the other sym- 
bols arc particle masses. We have already presented Fig. 
3, data at /3 = 5.5 on volume (12 x 8^) x 8. The other 
finite temperature data sets produce similar behavior. 
Fig. 7 shows spectra from (3 = 5.1, Fig. 8 shows spectra 
from f3 = 5.2, and Fig. 9 shows spectra from /3 = 6.0. 
The /3 = 5.1 data set shows a gap in the range of niq 
values which comes of the abrupt change in aniq with 
K (see Fig. 5). Here the points at large quark mass are 
confined, while in Figs. 8 and 9 all points are in the de- 
confined phase. In all these cases, the lightest Matsubara 
frequency is 7r/8 so the minimum meson mass according 
to Eq. (1) would be about 0.8. 

Scans of spectroscopy on 12^ lattices show similar be- 
havior to what we saw on the smaller lattices. Two exam- 
ples are shown in Fig. 10 and Fig. 11, f3 = 5.3 and 5.7. 
At P = 5.3 the three heaviest quark mass points are 
confined. At /? = 5.7 the two heaviest mass points are 
confined. The next point is on the transition and the 
rest arc deconfincd. These are not screening masses, 
but measurements performed in the temporal direction — 
ordinary spectroscopy. We have checked screening corre- 
lators at several of these points and they produce identi- 
cal results. 



C. Searches for curves of constant physics 

Motivated by experience with using lattice calculations 
for ordinary QCD phenomenology, we attempted to find 
lines of constant physics, along which we might attempt 
to make continuum predictions. We were not succcssftil 
in doing this. We describe briefly what we tried to do, 
and why it failed. 

A curve of constant physics is a manifold of points in 
the space of all bare couplings on which a dimensionless 
ratio of two dimensionful parameters (with the same mass 
dimension) is a constant. For our theory, this manifold is 
a line in the k, /? plane. Along that line, a lattice spacing 
can be defined from any dimensionful observable. (For 
example, we might find a line of constant mp/mv, and 
along that line, amy would give the lattice spacing, at 
the fixed ratio.) 

Our problem was twofold. First, we were squeezed 
between short-distance artifacts at small f3 and finite- 
volume effects at large /3. The former became apparent 
when the potential became essentially linear at all dis- 
tances, meaning that the Sommer scale had fallen to the 
neighborhood of the lattice spacing. The finite-vohune 
effects become apparent, of course, at the deconfinement 
transition in finite temperature and at the point where 
the linear potential disappears at finite spatial volume. 
As mentioned above, this left us with only a narrow strip 
in (/?, k) lying along and beneath the confinement tran- 
sition curve. 

The second problem was that we have two kinds of ob- 
servables, "gluonic" or "confinement" observables (such 

as To, cr, and Tc) and "mesonic" or "chiral" observables 
(such as mp, my and fp). Curves of constant physics 
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FIG. 8: Screening masses and fp for (3 = 5.2 on volume (12 x 8^) x 8. Labels as in Fig. 3. 



from confinement observables, such as lines of constant 
roy/a, simply did not look anything like curves from chi- 
ral physics, such as lines of constant mp/mv- Evidently 
there is no consistent way to define the lattice spacing in 
the region in question. 



V. SUMMARY AND FUTURE PROSPECTS 

The two main issues addressed in this work are the 

existence of separation between chiral and confinement 
scales and the manifestation of a supposed infrared fixed 



point in the massless theory. We have presented evidence 
against the former while the latter is supported by the 
behavior of the deconfinement transition curves as Nf is 
changed. 

Rather than seeing scale separation, we found only a 
single phase boundary in the k) plane for any given 
Nf. The weak-coupling side of this phase boundary shows 
parity doubling and a smooth connection to the region 
near Kc, indicating that there is no obstacle on the way to 
a chiral limit (or the Wilson-fermion equivalent) without 
spontaneous symmetry breaking. The smoothness of fp 
and its rapid decrease as Kc is approached is consistent 
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FIG. 10: Mass spectrum and fp for /3 = 5.3 on volume 12''. Labels as in Fig. 3. 



with this picture. 

We did not sec a chirally broken but deconfined phase. 
Our study does not, however, rule out such a phase near 
the would-be meeting point of the deconfinement and Kc 
curves. Here we were limited by the notorious numeri- 
cal difficulties associated with unimproved Wilson-clover 
fermions. Moreover, the lack of ehiral symmetry in the 
Wilson formulation (as opposed to staggered fermions 
[20, 22, 23]) leads to the absence of a true order parame- 
ter, such as the quark condensate, that could be extrap- 
olated to the massless limit. After all, all of the evidence 
we know for scale separation involves measuring the con- 
densate and the Polyakov loop, and seeing them order at 



different bare couplings. 

Now let us expand on the significance of the Kth curves 
for the question of the IRFP. Like all lattice gauge the- 
ories, our model possesses a strong coupling phase with 
a mass scale dictated by confinement. This region must 
be separated from the basin of attraction of an IRFP be- 
cause theories in this basin are conformal in the infrared, 
and thus cannot support a dynamically generated scale. 
This basin, in turn, has to be a subset of the aruq — 
(i.e., K = Kc) line, because a quark mass is a relevant 
coupling at a conformal fixed point. The basin includes 
f3 ~ oo but it might not extend all the way to /? = 0; it 
can end at a (UV-attractive) critical point at some (3* or 
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FIG. 11: Mass spectrum and fp for /3 

at a first order bulk transition. 

At finite volume — say, infinite spatial volume but 
nonzero Nf the system undergoes a confinement phase 
transition or crossover; the phase boundary starts at in- 
finite quark mass [k = 0) and extends into the diagram. 
Let us suppose that it extends to an intersection with 
the K = Kc hue, so that the massless theory exhibits 
the phase transition as well. (The alternative is that the 
phase boTindary curves off towards /? = at k < Kc.) 
At infinite mass our theory is a pure gauge theory so, as 
Nt increases, the transition will move to ever larger /3. 
At rUq = 0, however, the transition has to remain in the 
confinement region of the infinite volume theory, outside 
of the basin of the IRFP. This means that the deconfine- 
ment lines for different iV^'s have to have an accumulation 
point on the Kc line qualitatively resembling the situation 
shown in Fig. 1. This is the f3* critical point or the first 
order transition that bounds the IRFP's basin. 

A set of deconfinement lines that march up the Kc line, 
crossing the location of the putative IRFP, would have 
been inconsistent with the existence of an IRFP. Simi- 
larly, observing cliiral symmetry breaking or confinement 
at parameter values near the IRFP would have been an- 
other inconsistency. Our data point the other way. 

SU(2) gauge theory with Nf — 2 flavors of adjoint 
fermions is also a candidate for walking technicolor, and 
may instead possess an IRFP [16]. This theory shares 
the main features of our model. The SU(2) pure gauge 




0.0 0.2 0.4 0.6 

5.7 on volume 12'*. Labels as in Fig. 3. 

theory possesses a second order confinement transition 
at finite Nt and adjoint fermions preserve the associated 
Z2 center symmetry; hence the pattern of finite temper- 
ature transitions sweeping across the phase diagram at 
finite simulation volume is expected there, too. If the 
massless theory has an IRFP, then the zero-mass limit 
of the confinement transition must again terminate out- 
side of its basin of attraction, and the phase diagram will 
resemble our Fig. 1. 

We will return to the sextet theory with a study 
at smaller quark masses, which will become possible 
through the use of an improved lattice action. 
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